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Abstract 


We  have  undertaken  an  investigation  of  the  structure  and  dynamics  of  the  high-energy  ionic 
liquid,  l-hydroxyethyl-4-amino-l,2,4-triazolium  nitrate  (HEATN).  Both  experimental  and 
computational  methods  have  been  employed  to  understand  the  fundamental  properties,  characteristics, 
and  behavior  of  HEATN.  This  system  maintains  an  even  charge  separation,  as  assessed  by  both  the 
Mulliken  and  geodesic  derived  charges,  in  contrast  to  other  reports  on  triazolium  ionic  liquids.  The  MP2 
level  optimizations  find  six  dimer  and  five  tetramer  structures  and  show  the  significant  highly  hydrogen 
bonded  network  within  HEATN.  The  fragment  molecular  orbital  (FMO)  method  adequately  treats  this 
ionic  liquid  system  as  evidenced  by  the  small  relative  error  obtained  while  avoiding  the  prohibitive 
scaling  of  other  ab  initio  electronic  structure  methods.  The  solvation  dynamics  of  the  HEATN  system 
were  investigated  via  the  coumarin  153  probe  at  five  different  temperatures,  and  the  rotational  dynamics 
were  measured.  Comparisons  with  previously  studied  imidazolium  and  phosphonium  ionic  liquids  show 
surprising  similarity.  To  our  knowledge,  this  the  first  study  of  solvation  dynamics  in  a  triazolium  based 
ionic  liquid. 
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Introduction 


There  has  been  a  steeper  than  exponential  growth  in  the  number  of  publications  related  to  ionic 
liquids  since  the  year  2000. 1  The  interest  in  ionic  liquids  is  undoubtedly  owing  to  the  environmentally 
friendly  nature  of  these  compounds.  '  Ionic  liquids  can  be  liquid  at  room  temperature  and  are  molten 
salts  commonly  composed  of  an  organic  cation  and  an  inorganic  anion.  These  low-melting  salts  are 
considered  green  solvents,  in  contrast  to  volatile  organic  compounds  (VOCs),  due  to  their  high  chemical 
and  thermal  stabilities,  negligible  vapor  pressure,5  and  high  recoverability  and  reusability.67  Ionic 
liquids  can  be  customized  by  varying  the  identity  of  the  substituents  on  the  cation  or  adapting  the 

o 

moiety  of  the  anion.  This  versatility  and  these  unique  properties  have  piqued  the  interest  of  the 
scientific  community,  and  a  variety  of  fundamental  studies9"15  have  been  performed  on  them  to  obtain  a 
better  understanding  of  their  characteristics. 

One  interesting  application  of  ionic  liquids  is  in  the  high  energy  density  materials  (HEDM) 
community  where  many  different  molten  salts  have  been  explored  and  reported.16  It  is  important  to 
understand  what  characteristics  make  ionic  liquids  suited  to  HEDM  applications.  This  very  problem  was 
explored  computationally  in  a  series  of  studies,  "  as  well  as  experimentally  in  a  few  reports.  These 
studies  show  that  triazolium  ionic  liquids  appear  to  have  more  suitable  characteristics  for  HEDM 
application  than  other  nitrogen  rich  ionic  liquid  species. 

A  recent  report  by  Drake  and  co-workers  on  low-melting,  potentially  energetic,  salts  lists  the 
desired  characteristics  for  these  salts,  as  well  as  potential  benefits  that  might  be  derived  from  using  them 
in  energetic  materials  applications.  The  same  authors  have  discussed  the  l-R-4-amino-l,2,4-triazolium 
family  of  salts"  that  are  of  interest  owing  to  their  high  nitrogen  content.  One  interesting  room 
temperature  ionic  liquid  in  this  family  is  l-hydroxyethyl-4-amino-l,2,4-triazolium  nitrate  (HEATN), 
depicted  in  Figure  1.  To  use  this  ionic  liquid  in  various  applications  it  is  beneficial  to  employ  both 
experimental  and  computational  methods  to  understand  its  fundamental  properties,  characteristics,  and 
behavior. 
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This  present  work  focuses  on  the  physical  chemistry  of  the  high-energy  ionic  liquid  system 
HEATN.  The  HEATN  system  was  recently  investigated  using  molecular  dynamics  simulations  for 
prediction  of  bulk  properties.”  The  studies  presented  here  assess  this  ionic  liquid  at  a  molecular  level, 
by  using  ab  initio  electronic  structure  methods,  to  gain  insight  into  the  structure,  energetics  and  charge 
delocalization  in  the  component  ions.  These  properties  have  been  shown  to  be  of  significance  in  relation 
to  the  physical  and  chemical  properties  of  ionic  liquids,  and  a b  inito  electronic  structure  methods  have 
been  shown  to  give  excellent  agreement  with  X-ray  structures  of  triazole  salts  giving  bond  lengths 
within  0.003  nm  of  the  experimental  values  in  the  case  of  the  common  3, 4, 5-triamino- 1,2, 3-triazole 
salts.  Performing  electronic  structure  calculations  on  large  systems  of  ionic  liquids  can  be  difficult 
given  the  fact  that  ab  initio  electronic  structure  theory  methods  scale  ~N4  or  worse,  depending  on  the 
chosen  method  and  if  electron  correlation  is  included,  where  N  measures  the  size  of  the  system.  This 
means  that  correlated  electronic  structure  methods  become  computationally  prohibitive  for  very  large 
systems  with  more  than  100  heavy  atoms.  A  typical  approach  for  large  systems  is  to  treat  sections  of 
higher  interest  using  quantum  mechanics  (QM)  and  other  section  using  molecular  mechanics  (MM). 
These  combined  QM/MM  methods  have  had  some  success  for  large  systems  yet  at  times  important 
contributions  from  parts  of  the  system  may  be  treated  by  classical  model  potentials  having  a  significant 
impact  on  the  QM  region  leading  to  the  sacrificing  of  high  accuracy  to  save  time  and  resources. 
Alternatively,  breaking  the  large  system  up  into  smaller  sections  and  treating  these  fragments  with 
particular  levels  of  electronic  structure  theory  has  the  advantage  of  being  virtually  fully  quantum 
mechanical.  Fragmentation  methods  have  a  number  of  distinct  advantages  without  loosing  accuracy  for 
larger  systems.”  ’  The  fragment  molecular  orbital  (FMO)  method”  "  splits  the  system  up  into 
fragments  and  computes  each  fragment  (monomer)  in  the  electrostatic  potential  of  the  other  fragments. 
The  FMO  method  is  suited  for  larger  systems  and  has  been  applied  to  systems  such  as  polypeptides  and 
proteins,  and  even  a  heavy  metal  system  containing  3596  atoms.  Farge  systems  can  be  treated  at  a 
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high  level  of  theory  allowing  valuable  insight  into  large  systems  of  ionic  liquids.  The  use  of  the  FMO 
method  with  ionic  liquid  systems  has  been  limited  and  is  explored  in  this  report. 

Additional  studies  used  to  investigate  the  HEATN  system  and  its  dynamics  include  the  use  of 
fluorescence  spectroscopy  to  compare  this  high-energy  ionic  liquid  to  other  more  commonly  used  ionic 
liquids.  Steady-state  and  time-resolved  fluorescence  measurements  were  performed  to  assess  how 
HEATN  behaves  as  a  solvent  at  short  time  scales.  In  particular,  solvation  dynamics  experiments  were 
carried  out  using  the  fluorescent  probe  coumarin  153  (C153)  in  HEATN.  The  structure  of  C153  is  given 
in  Figure  1.  C153  is  an  ideal  probe,  and  has  been  used  extensively  in  a  large  number  of  solvation 
dynamics  studies  in  various  media.34  50  This  permits  facile  comparison  of  the  HEATN  ionic  liquid  with 
other  systems.  A  summary  of  the  various  solvation  dynamics  studies  in  ionic  liquids  has  recently  been 

c  i  co 

published  as  well  a  report  with  some  of  the  applications  of  such  explorations.  Up  till  now  there  have 
been  no  explorations  of  the  solvation  dynamics  of  triazolium  ionic  liquid  systems.  The  solvation 
dynamics  of  the  HEATN  system  were  investigated  via  the  C153  probe  at  five  different  temperatures, 
and  the  rotational  dynamics  were  measured  by  exploiting  the  anisotropic  nature  of  the  probe’s 
fluorescence  through  this  viscous  ionic  liquid  and  comparisons  with  previously  studied  imidazolium  and 
phosphonium  ionic  liquids  are  presented. 

Materials  and  Methods 
Computational  methods. 

Molecular  structures  were  obtained  by  performing  geometry  optimization  calculations  using 

co 

second-order  Moller-Plesset  perturbation  theory  (MP2)  with  a  Hartree-Fock  reference.  All  calculations 
were  made  using  the  6-31++G(d,p)  Pople  type  basis  set.  Hessians  (matrices  of  the  second  derivatives  of 
the  energies)  were  determined  to  ensure  each  stationary  point  was  truly  a  minimum  or  a  transition  state. 
Population  analysis  was  completed  using  Mulliken’s  method54'57  and  geodesic  derived  charges.58 

oo 

Studies  of  three  cation-anion  pairs  (hexamers)  were  carried  out  using  the  FMO  method.”  The  fragment 

molecular  orbital  method  provides  a  way  of  treating  large  systems  with  high  accuracy  by  overcoming 
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the  prohibitive  scaling  associated  with  ab  initio  electronic  structure  theory  methods.  This  is 
accomplished  by  dividing  the  large  molecular  system  of  interest  into  fragments  (each  called  monomers) 
and  executing  the  now  feasible  ab  initio  calculations  on  the  fragments  and  their  interaction  with  another 
fragment  (dimers)  or  two  other  fragments  (trimers)  in  the  presence  of  the  other  monomers. “  Such  that 

N  N 

the  total  energy  of  a  fragmented  dimer  system  is  given  as  follows:  E  =  £e, 

l  I>J 

where  the  energy  (E)  is  given  by  the  standard  SCF  method  for  each  monomer  (7)  for  a  total  of  N 
fragments  in  the  electrostatic  field  of  the  additional  N-l  fragments,  and  dimer  (77)  calculated  in  the 
electrostatic  field  of  the  other  N-2  fragments.  The  FMO  method  has  been  enhanced  by  the  use  of  the 
generalized  distributed  data  interface  (GDDI)  allowing  two  levels  of  parallelization,59  and  is  available  in 
a  multilevel  format  allowing  one  to  assign  fragments  to  particular  levels  so  that  you  can  treat  these 
different  levels  with  different  electronic  structure  methods  and/or  basis  sets.60  With  this  ionic  liquid 
system  one  molecular  ion  was  assigned  to  a  fragment,  thereby  removing  the  need  to  fragment  bonds  as 
is  common  with  large  covalently  bonded  systems.  For  covalently  bonded  systems  there  are  some 

no 

guidelines  that  can  help  to  obtain  the  best  level  of  accuracy. “  The  FMO  calculations  were  carried  out 
using  the  two-body  correction  in  conjunction  with  MP2  (FM02-MP2)  as  well  as  the  three-body 
correction  with  MP2  (FM03-MP2).  All  calculations  were  performed  using  GAMESS,61'62  and  were 
visualized  when  possible  with  MacMolPlt.63 

Materials. 

Coumarin  153  (C153)  (Exciton  Inc.,  Dayton,  OH)  was  used  as  received.  4-amino- 1,2,4-triazole 
was  purchased  from  Sigma.  The  production  of  l-substituted-4-amino-l,2,4-triazolium  nitrate  ionic 
liquid  molecules  through  a  metathesis  of  heterocyclic  halide  and  nitrate  salts  is  well  established.22’64 
The  specific  synthesis  of  l-hydroxyethyl-4-amino-l,2,4-triazolium  bromide  (HEATB)  and  1- 
hydroxyethyl-4-amino-l,2,4-triazolium  nitrate  (HEATN)  from  metathesis  of  HEATB  with  silver  nitrate 

follows  the  procedure  described  by  Drake,  Hawkins  and  Tollison.65 
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The  synthesis  of  the  other  ionic  liquids  shown  in  Figure  1  is  described  elsewhere  10  and  were 
decolorized  according  to  the  procedure  described  elsewhere.66  Briefly,  the  impure  ionic  liquid  was 
diluted  and  was  allowed  to  elute  through  a  column  packed  with  celite,  silica  gel  and  activated  charcoal. 
After  the  excess  solvent  was  removed,  the  ionic  liquid  was  dried  under  vacuum  with  mild  heating.  The 
water  content  of  the  HEATN  sample  was  determined  using  a  Mettler  Toledo  DL  39  coulometric  Kail 
Fischer  titrator.  The  viscosity  of  the  HEATN  sample  was  measured  at  each  of  the  five  temperatures  ± 
0.1  °C  using  a  ViscoLab  4000  viscometer  from  Cambridge  Applied  Systems. 

Steady-state  measurements. 

Steady-state  absorption  spectra  were  obtained  using  a  Hewlett-Packard  8453  UV-visible 
spectrophotometer  with  1-nm  resolution.  Steady-state  fluorescence  measurements  were  taken  using  a 
SPEX  flurormax-4  spectrofluorometer  (HORIBA  Jobin  Yvon)  with  1-nm  resolution.  The  emission 
spectra  were  corrected  for  detector  response  and  the  lamp  spectral  intensity.  All  emission  spectra  were 
obtained  with  an  excitation  wavelength  of  407  nm  and  a  3  nm  band-pass.  A  3  mm  path  length  quartz 
cuvette  was  used  for  all  measurements. 

Time-resolved  measurements. 

Fluorescence  lifetime  measurements  were  made  using  the  time-correlated  single-photon 
counting  (TCSPC)  apparatus  described  previously9’50  a  summary  is  included  here.  Using  a  home-made 
mode-locked  Ti:sapphire  oscillator  pumped  at  532  nm  by  a  Nd:VC>4  Millennia  (Spectra-Physics)  laser 
femtosecond  pulses  were  produced  which  were  tunable  from  780  nm  -  900  nm  having  a  repetition  rate 
of  82  MHz.  With  the  output  selected  at  814  nm  from  the  Ti:Sapphire  oscillator,  the  repetition  rate  was 
reduced  to  8.8  MHz  via  the  use  of  a  Pockels  cell  (model  350-160  from  Conoptics  Inc.)  and  was 
subsequently  frequency  doubled  through  the  use  of  a  harmonic  generator  (model  TP-2000B  from  U- 
Oplaz  Technologies).  The  produced  blue  light  having  a  wavelength  of  407  nm  was  used  as  the  primary 
excitation  source  for  all  time-resolved  studies.  After  the  harmonic  generator  a  half-wave  plate  was 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


placed  in  front  of  a  vertical  polarizer  to  ensure  the  polarization  of  the  excitation  light.  The  fluorescence 
was  collected  at  a  90°  angle  to  the  excitation  source  and  then  passed  through  an  emission  polarizer  set  to 
the  magic  angle  (54.7°)  with  respect  to  the  vertical  excitation  light.  In  the  case  of  the  anisotropy  studies 
this  polarizer  was  set  to  0°  and  90°  with  respect  to  the  vertical  excitation  light  respectively.  A  425  nm 
cut-off  filter  was  placed  before  a  multichannel  plate  (MCP)  (Hamanatsu)  to  help  remove  any  unwanted 
scattered  excitation  light.  When  performing  the  solvation  dynamics  studies,  a  monochromator  was 
placed  before  the  MCP  to  ensure  the  collection  at  each  selected  wavelength.  The  signal  from  the  MCP 
detector  was  amplified  and  sent  to  a  Becker  &  Hickl  photon  counting  module  (model  SPC-630).  The 
instrument  response  function  had  a  full  width  at  half  maximum  (FWHM)  of  ~45  -  50  ps.  The  parallel 
and  perpendicular-polarized  fluorescence  anisotropy  decay  curves  were  collected  in  a  22  ns  time- 
window  and  fitted  simultaneously  following  the  method  described  by  Cross  and  Fleming.67  This 
method  takes  full  advantage  of  the  statistical  properties  of  the  measured  curves.  The  solvation  dynamics 
studies  involved  recording  fifteen  single  wavelength  decays  from  490  nm  to  630  nm  in  10  nm  intervals, 
in  a  time-window  of  9  ns.  Typically,  1024  channels  of  data  were  collected  with  65,000  counts  in  the 
peak  channel.  The  wavelength-resolved  fluorescence  transients  were  fit  to  sums  of  exponentials 
(typically  2  or  3,  as  necessary  to  fit  the  data),  and  time-resolved  emission  spectra  (TRES)  were 
reconstructed  as  described  in  previous  reports.9'50 

The  traditional  approach  of  fitting  the  time -resolved  emission  spectra  to  a  log-normal  function, 
9.34,68  from  we  extract  the  peak  frequency  v(t)  as  a  function  of  time  was  used.  The  solvation 


dynamics  is  described  by  the  following  normalized  correlation  function: 

C(t)_  v(t) -v(co) 
v("0")  —  v(co) 


(1) 


Because  C(t)  is  a  normalized  function,  the  accurate  determination  of  C(l)  depends  upon  accurate  values 

for  v(“0”)  and  v(oo).  v(“0”)  is  the  frequency  at  zero-time,  estimated  using  the  method  of  Fee  and 

Maroncelli,69  who  have  described  a  robust,  model  independent,  and  simple  procedure  for  generating  this 
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“zero-time”  spectrum,  v(“0”)  which  represents  the  emission  spectrum  expected  prior  to  any  solvent 
relaxation  but  after  complete  intramolecular  vibrational  redistribution.  The  validity  has  been  checked 
using  a  different  method  for  estimating  the  “zero-time”  reorganization  energy  in  a  previous  report.10 
v(oo)  is  (usually70,71)  the  frequency  at  infinite  time,  obtained  from  the  maximum  of  the  steady  state 
spectrum.  (This  is  not,  however,  true  in  the  case  of  very  slowly  relaxing  solvents,  as  has  been 
demonstrated  in  the  case  of  certain  ionic  liquids  ’  ’  :  where  the  emission  spectrum  at  ~3  times  the 
fluorescence  lifetime  of  the  probe  is  reel-shifted  to  that  of  the  equilibrium  spectrum.)  The  v(t)s  are 
determined  from  the  maxima  of  the  log-normal  fits  of  the  time  resolved  emission  spectra.  In  most  of  the 
cases,  however,  the  spectra  are  broad,  so  there  is  some  uncertainty  in  the  exact  position  of  the  emission 
maxima.  Thus,  the  range  of  the  raw  data  points  in  the  neighborhood  of  the  maximum  to  estimate  an 
error  for  the  maximum  obtained  from  the  log-normal  fit  was  considered.  Depending  on  the  width  of  the 
spectrum  (i.e.,  “zero-time”,  steady-state,  or  time-resolved  emission  spectrum),  it  was  determined  that  the 
typical  uncertainties  are  as  follows:  “zero-time”  ~  steady-state  (~  ±  100  cm'1)  <  time-resolved  emission 
(~  ±  200  cm  ').  These  uncertainties  were  used  compute  error  bars  for  the  C(l)  graph  (giving  maximum 
error  of  ±0.04).  Finally,  in  generating  the  C(t )  curve,  the  first  point  was  obtained  from  the  “zero-time” 
spectrum.  The  second  point  was  taken  at  the  maximum  of  the  instrument  response  function.  The 
fractional  solvation  was  also  computed  at  50  ps  by  taking  f5ops=l-C(t=50ps).  The  experiments  were  also 
carried  out  at  five  different  temperatures  (25  °C,  35  °C,  45  °C,  55  °C  and  65  °C)  that  were  obtained 
using  a  variable  temperature  sample  cell  holder. 

Results  and  Discussion 

HEATN  Structures.  A  previous  ab  initio  investigation17  of  triazolium  ionic  liquid  systems 

showed  a  tendency  for  proton  transfer  between  the  cation  and  the  anion  to  form  neutral  pairs.  This 

tendency  was  investigated  within  the  HEATN  system  by  performing  geometry  optimizations  on  the 

different  pairs  of  cations  and  anions.  The  geometric  structures  of  both  single  cation-anion  pairs  (dimers) 

and  two  cation-anion  pairs  (tetramers)  have  been  investigated.  Six  structures  labeled  I,  II,  III,  IV,  V,  and 
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VI  were  found  for  the  HEATN  dimer  (Table  1)  at  the  MP2  level.  They  are  ordered  with  respect  to  their 
relative  energies  with  structure  “I”  having  the  lowest  energy.  The  relative  energies  can  be  seen  in  Table 
1  along  with  the  hydrogen  bond  distances,  Mulliken  charges,  and  geodesic  derived  charges.  The  charge 
delocalization  of  each  of  these  systems  shows  a  very  even  charge  separation  between  the  cation  and  the 
anion  as  expected.  Typically  the  geodesic  charges  are  considered  to  be  superior  to  the  Mulliken  charges; 
yet  good  agreement  between  these  two  methods  is  seen  in  most  of  these  cases.  Each  of  these  structures 
show  two  hydrogen  bonds  between  the  cation  and  anion  giving  evidence  that  the  HEATN  system  is  a 
highly  hydrogen  bonded  system.  This  is  in  part  by  design,  as  this  liquid  was  chosen  due  to  its  potential 
as  a  perspective  fuel.  Having  a  higher  density,  which  the  hydrogen  bonds  may  contribute  to,  is  one  of 
the  desirable  properties  of  high  energy  materials  as  this  allows  more  energy  per  cubic  centimeter.  Since 
no  geometries  were  found  where  the  anion  is  to  the  side  of  the  cation  it  seems  that  the  hydroxyethyl  side 
chain  on  the  cation  allows  for  the  carousal  of  the  nitrate  anion. 

In  addition,  the  HEATN  system  was  investigated  at  the  tetramer  (four  molecules)  level  using  the 
MP2  method.  The  geometries  obtained  can  be  seen  in  Figure  2.  The  relative  energies,  hydrogen  bond 
distances,  and  both  the  Mulliken  and  geodesic  charges  are  given  in  Table  2.  Interestingly,  these 
structures  show  less  sharing  of  the  anions  between  the  cations.  The  cations  seem  to  prefer  to  interact 
with  only  one  anion  and  are  correlated  via  hydrogen  bonds  to  only  a  single  anion.  The  charge  separation 
is  very  even  in  the  lowest  energy  structure  found  while  the  other  structures  show  a  less  even 
distribution. 

The  FMO  method  was  used  to  evaluate  its  feasibility  for  use  with  larger  (hexamer  and  higher) 

charged  ionic  systems  such  as  this  ionic  liquid.  These  structures  serve  as  a  great  reference  point  to 

explore  how  both  the  FM02-MP2  and  the  FM03-MP2  method  handle  this  highly  energetic  ionic  liquid 

system.  Applying  these  FMO  methods  to  the  MP2  geometries  yields  the  results  as  seen  in  Table  3.  As 

expected,  the  three-body  correction  (FM03-MP2)  does  a  better  job  of  recovering  the  MP2  results. 

Comparing  the  results  from  the  FM02  and  FM03  methods  it  is  immediately  seen  that  they  both  recover 
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the  MP2  energy  well  with  the  largest  absolute  difference  being  <  1  kcal/mol  in  the  FM02  case  and  < 
0.30  kcal/mol  using  the  FM03  method.  While  in  general  better  results  are  obtained  using  the  three -body 
correction  it  is  often  not  employed  with  smaller  systems  because  of  the  additional  computer  time  and 
resources  required.  In  this  particular  case  the  MP2  method  used  71.55  MB  of  RAM  while  the  FM02- 
MP2  and  FM03-MP2  methods  used  5.07  MB  and  11.47  MB  respectively.  It  seems  that  as  the  system 
grows  the  extra  cost  of  using  the  FM03-MP2  method  will  be  outweighed  by  the  more  accurate  results 
obtained.26  The  charges  found  using  both  the  Mulliken  and  geodesic  charge  schemes  derived  from  the 
FM02  and  FM03  single  point  calculations  show  good  agreement  to  those  found  as  derived  from  the 
MP2  method.  The  geodesic  charges  show  that  there  is  an  acceptable  treatment  of  the  charges  by  the 
Mulliken  scheme,  at  least  in  the  case  of  this  ionic  liquid. 

The  accuracy  using  the  FMO  methods  can  be  increased  by  making  a  few  different  changes. 
Larger  systems  can  be  more  accurately  (in  comparison  to  the  MP2  results)  calculated  by  including  more 
molecules  per  fragment.  This  would  increase  the  costs  of  such  calculations  but  very  large  systems  can 
achieve  additional  accuracy  this  way.  Fragment  methods  must  normally  be  tailored  to  specific 
applications  to  achieve  higher  accuracy.  That  is,  particular  schemes  of  how  many  molecules  per 
fragment  and  how  to  divide  molecules  into  fragments  may  vary  for  different  systems.  More  work  needs 
to  be  done  on  larger  systems  of  HEATN  molecules  to  determine  how  exactly  the  HEATN  system  should 
be  partitioned  into  fragments  to  achieve  the  most  reasonable  level  of  accuracy  for  large  systems  of 
HEATN. 

Solvation  Dynamics  of  C153  in  HEATN 

HEATN  was  investigated  using  the  fluorescent  probe  molecule  C153.  The  maximum  of  the 

emission  spectrum  maximum  is  at  546  nm  and  is  typical  of  C153  in  other  solvents.  There  was  no 

temperature  dependent  shift  of  the  absorption  and  emission  spectra  of  C153  in  HEATN.  In  order  to 

construct  the  solvation  correlation  function  C(t),  time-resolved  fluorescence  decays  were  collected  at 

different  wavelengths  spanning  the  entire  range  of  the  emission  spectrum  of  the  probe  as  shown  in 
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Figure  3.  The  decays  on  the  blue  end  were  characterized  by  a  faster  response  whereas  those  collected 
at  longer  wavelengths  were  accompanied  by  a  rise  time  which  was  due  to  the  population  flux  from 
initially  excited  state  to  relaxed  state.  This  trend  of  wavelength  dependent  decays  is  a  signature  of  local 
solvation  dynamics.  Time -resolved  emission  spectra  at  50  ps  and  7.05  ns  of  C153  in  HEATN  are 
presented  in  Figure  4  at  two  representative  temperatures.  The  corresponding  “steady-state”  and  “zero¬ 
time”  spectra  are  also  included.  Note  that  unlike  at  65°C,  the  7.05  ns  spectrum  crosses  the  steady-state 
spectrum  at  25°C,  indicating  that  this  spectrum  does  not  adequately  report  on  the  completed  solvation  of 
C153,  as  would  usually  be  expected  for  highly  viscous  solvents.  ’  ’  The  spectrum  recorded  at  7.05  ns 
was  red-shifted  by  -270  cm  1  relative  to  the  steady-state  spectrum,  demonstrating  that  solvation  is 
slightly  slower  than  the  population  decay  of  the  Si  state  of  Cl 53.  This  phenomenon  has  been  previously 
observed  in  the  case  of  the  slowly  relaxing  phosphonium  ionic  liquids.  In  previous  studies,  it  was 
found  that  the  amount  of  rapid  solvation  and  the  average  solvation  time  decrease  and  increase, 
respectively,  as  the  time-window  of  the  experiment  is  increased.  This  phenomenon  of  solvation  slower 
than  the  population  decay  of  the  Si  state  has  also  been  noticed  elsewhere  for  related  systems.  ’  In 
order  to  compare  the  solvation  in  the  ionic  liquid  media,  we  attempted  to  study  the  same  in  parent 
compound  which  is  the  neutral  4-amino- 1,2,4-triazole  (melting  point  -  85  °C).  A  solvent  induced 
spectral  shift  was  not  observed  in  the  molten  state  at  85  °C  with  the  TCSPC  apparatus,  suggesting  a  very 
fast  response  for  the  triazole,  which  is  consistent  with  its  low  viscosity  (~  6.2  cP)  at  the  melting 
temperature. 

The  spectral  shift  dynamics  are  represented  by  the  decay  of  the  solvation  correlation  function 
C(t)  in  Figure  5.  The  C(t) s  are  constructed  at  five  temperatures  (25,  35,  45,  55  and  65°C)  and  the 
solvation  parameters  are  collated  in  Table  4.  At  all  temperatures,  the  solvation  dynamics  were  well 
described  by  biphasic  decays,  with  an  initial  faster  component  preceded  by  a  slower  response.  At  25°C, 
almost  >30%  of  the  solvation  is  completed  within  instrumental  time  resolution  (~50ps),  which  is 
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determined  from  the  fractional  solvation,  /sops  listed  in  Table  4.  This  accounts  for  the  solvent  response 
missed  due  to  the  finite  time -resolution  of  the  TCSPC  apparatus.  The  missing  Stokes-shift  was 
accounted  for  by  using  the  “zero-time”  spectra  estimated  by  the  method  of  Fee  and  Maroncelli.69 
Recently  Kerr-gated  emission  (KGE)  ’  and  fluorescence  upconversion  techniques  '  were  used  to 
measure  spectral  dynamics  in  ionic  liquids  with  must  shorter  time -resolution  ranging  from  200  fs  -  500 
fs.  Both  the  Maroncelli  and  Vauthey  groups  have  reported  a  small  but  significant  amount  of 
subpicosecond  component  in  the  case  of  solvation  in  ionic  liquids.  Maroncelli  and  coworkers  used  a 
combination  of  KGE  and  TCSPC  data  to  resolve  complete  solvation  response  in  ionic  liquids  using  4- 
dimethylamino-4’-cyanostilbene  (DCS)  as  a  probe.  '  The  observed  response  functions  were  well 
described  by  biphasic  decays,  consisting  of  both  a  sub-picosecond  component  of  -10-20%  amplitude 
and  a  dominant  slower  component  relaxing  over  a  few  picoseconds  to  several  nanoseconds.  The  faster 
component  was  correlated  to  associated  inertial  characteristics  of  the  constituent  ions,  and  the  slower 
component  to  solvent  viscosity.  The  missing  portion  of  the  initial  relaxation  contains  both 
subpicosecond  inertial  contributions  as  well  as  components  relaxing  on  the  1-10  ps  time  scale  '  and  this 
interpretation  was  also  corroborated  by  simulations  performed  by  Kobrak.76 

As  expected  the  solvation  times  decreased  with  increasing  temperature.  This  was  most  likely 

because  of  the  decrease  in  the  viscosity  of  the  ionic  liquid  allowing  the  neighboring  solvent  molecules  to 

reorganize  faster  around  the  excited  state  dipole  of  C153.  Figure  6a  shows  the  dependence  of  solvation 

time  on  viscosity  in  different  types  of  ionic  liquids  namely  those  based  on  imidazolium,  phosphonium 

and  triazolium  cations.  The  solvation  times  of  C153  in  HEATN  show  a  reasonable  correlation  with 

viscosity,  whereas  in  imidazolium  based  ionic  liquids  the  correlation  is  not  as  good;  and  marked 

deviations  are  observed  for  the  phosphonium  ionic  liquids.  This  observation  is  consistent  with  those 

reported  by  Maroncelli  and  coworkers,  ’  where  they  attributed  the  deviations  in  phosphonium  ionic 

liquids  to  the  enormous  size  of  the  cation  moiety.  Solvation  in  phosphonium  ionic  liquids  is  much 

slower  than  that  of  the  imidazolium  and  HEATN  systems.  For  example,  although  (Ce)3Ci4P+  Br"  ‘ 77  is 
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isoviscous  (-270  cP)  with  HEATN  (at  35°C),  the  solvation  is  5  times  slower  in  the  former.  This 


observation  can  be  interpreted  in  terms  of  the  van  der  Waal  volumes  of  the  constituent  ions.  The  radii 
of  the  cations  and  anions  as  listed  in  Table  4  were  estimated  using  Edward’s78  and  Bondi’s79  van  der 

o 

Waal  increment  methods.  The  radii  decrease  in  the  order  of  phosphonium  (-5. 3 A)  >  imidazolium 

o  o 

(-3.4A)  >  triazolium  (-2.9A).  This  probably  also  explains  the  2-fold  faster  average  solvation  times  in 
HEATN  compared  to  the  isoviscous  imidazolium  based  ionic  liquids.  It  is  noteworthy  in  this  context 
that  in  all  the  imidazolium  based  ionic  liquids  studied  so  far,  the  completely  relaxed  state  or  equilibrium 
condition  is  represented  adequately  by  the  steady-state  spectrum  of  the  probe,  however  this  was  not  the 
case  in  HEATN  under  ambient  temperature. 

As  discussed  above,  besides  viscosity  (q),  the  size  of  the  individual  ions  strongly  influences  the 
solvent  response.  In  order  to  investigate  how  the  ion  size  might  affect  the  observed  solvation  times,  we 
have  performed  different  analyses  motivated  by  the  methods  of  Maroncelli  and  coworkers,  which  were 
used  to  establish  correlations  of  the  logarithm  of  the  solvation  times  (ts)  to  various  measured  and 
calculated  properties.  A  plot  of  log  <is>  versus  q/T  did  not  provide  substantial  improvement  compared 
to  a  plot  of  log  <xs>  against  q  only.  Consistent  with  the  reports  of  Jin  et  al.,  77  a  much  better  correlation 
was  found,  as  shown  in  Figure  6b,  in  which  the  solvation  time  varies  as  (r//T)A  R+  ,  where  R+  is  the 
radius  of  the  cation,  and  A  and  B  are  constants  determined  to  be  0.5  and  4.5  by  multiple  regression 
analysis.  This  correlation  was  slightly  improved  when  the  anion  size  ( R  .)  is  incorporated  as  shown  in 
Figure  6c.  Thus,  the  dependence  of  solvation  time  could  be  best  described  as  follows: 

<  >oc  <2) 

The  diffusion  times  of  the  cation  and  anion  of  HEATN,  which  are  the  times  to  move  a  root  mean 
square  distance  equal  to  its  radius,  were  calculated  using  the  relation  t±  =  ///?;'  /  f±kBT ,  where  is  the 
Boltzman  constant  and  /  is  0.53  and  0.74  for  ‘+’  and  which  refer  to  the  cation  and  anion 
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respectively.  The  diffusion  times  are  presented  in  Table  4.  The  correlation  of  solvation  time  with  the 
diffusion  time  of  the  cation  and  anion,  as  shown  in  Figure  7,  indicates  that  the  contribution  of  the 
reorientational  motion  of  ions  to  the  overall  relaxation  process  cannot  be  neglected. 

Rotational  Dynamics  of  Cl 53  in  HEATN 

The  fluorescence  anisotropy  of  Cl 53  in  HEATN  was  measured  at  five  different  temperatures 
and  decreased  with  temperature.  Unlike  in  normal  polar  solvents,  the  anisotropy  decay  was  non¬ 
exponential  at  all  temperatures.  The  limiting  anisotropy  (ro)  and  rotational  times  are  collated  in  Table  4. 
Figure  8  gives  the  variation  of  rotational  time  (rrot)  of  Cl 53  in  HEATN  as  a  function  of  temperature  and 
viscosity  in  the  form  of  p/T.  This  plot  shows  a  good  correlation  in  the  form<  Trot  >cc  A(rj /T)o  s ,  where 

A  is  a  constant  determined  from  the  linear  fit  to  be  1.1.  Ionic  liquids  have  been  reported  to  show  similar 

11 

correlation  as  established  by  normal  polar  solvents. 

The  experimentally  measured  anisotropy  decay  times  of  Cl 53  in  HEATN  have  been  analyzed 
within  the  framework  of  Stokes-Einstein-Debye  (SED)  hydrodynamic  theory.  According  to  this  theory, 
the  rotational  diffusion  of  a  medium  sized  solute  molecule  in  a  solvent  continuum  is  assumed  to  occur 
by  small  step  diffusion  and  its  reorientation  time  and  is  related  to  the  macroscopic  viscosity  of  the 

o  i  oo 

solvent  by  the  following  relation  ’  “ 

(3) 

Kb1 

The  dashed  lines  in  Figure  8  are  the  predictions  of  hydrodynamic  models,  made  by  assuming  C153  is  of 
ellipsoidal  shape  with  semi-axis  dimensions  of  2.0,  4.8,  and  6.1  A  and  the  relation  where  V  is  the 

°  o 

solute  van  der  Waals  volume  (246  A  ),  /is  a  shape  factor  (1.71),  and  C  is  a  “rotational  coupling  factor” 
which  signifies  the  extent  of  coupling  between  the  solute  and  the  solvent.  For  stick  boundary  conditions 

oo 

C  =1,  and  for  slip  boundary  conditions  C=0.24.  All  of  the  rotation  times  in  HEATN  measured  at 

different  temperatures  fall  between  these  two  limiting  predictions.  Deviations  from  the  hydrodynamic 

predictions  can  be  estimated  from  the  value  of  the  rotational  coupling  factor  (Crot)  using  the  relation 
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(4) 


where  zstick  is  the  stick  hydrodynamic  prediction.  The  values  of  Crot  are  listed  in  Table  4.  The  average 

77 

value  in  HEATN  is  0.55,  whereas  those  in  conventional  polar  solvents  were  reported  to  be  0.57+0.09. 
Thus  the  coupling  of  solute  and  solvent  between  Cl 53  and  HEATN  is  reasonably  comparable  to 
conventional  polar  solvents.  These  results  are  consistent  with  those  of  Maroncelli  and  coworkers  who 
reported  that  all  ionic  liquids,  besides  bulky  ones  based  on  phosphonium,  fall  within  the  scatter  of  the 
polar  solvents.  Castner  and  coworkers  have  recently  studied  solvation  and  rotational  dynamics  using 
Cl 53  over  a  range  of  temperatures  from  5  to  80°C,  in  ionic  liquids  based  on  ammonium  and 
pyrrolidinium  cations.  They  reported  small  amplitudes  of  very  long  orientation  relaxation  time 
constants  on  the  order  of  100  ns,  which  were  missing  in  the  experiments  reported  on  herein.  All  of  the 
Cl 53  rotational  time  constants  in  HEATN  ranged  from  2  -  20  ns  and  are  in  agreement  with  the  reports 
by  Maroncelli  and  coworkers.  Observation  of  dynamics  in  the  time-scale  of  100  ns  may  be  unreliable 
where  the  lifetime  of  C153  is  ~3-5  ns. 


Conclusions 

The  structure  and  dynamics  of  the  high-energy  ionic  liquid  system,  l-hydroxyethyl-4-amino- 
1,2,4-triazolium  nitrate  (HEATN),  were  investigated.  This  system  maintains  an  even  charge  separation 
assessed  by  both  the  Mulliken  and  geodesic  derived  charges.  There  were  no  instances  where  proton 
transfer  was  found  in  contrast  to  other  reports  on  triazolium  ionic  liquids.  The  MP2  level  optimizations 
find  six  dimer  and  five  tetramer  structures  and  show  the  significant  highly  hydrogen  bonded  network 
within  HEATN.  The  FMO  method  adequately  treats  this  ionic  liquid  system  as  evidenced  by  the  small 
relative  error  obtained.  Further  calculations  using  the  FMO  method  on  this  ionic  liquid  system  are 
expected  to  give  similar  accuracy  while  using  less  computer  resources  and  avoiding  the  prohibitive 
scaling  of  other  ab  initio  methods. 

Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Besides  the  structural  characterization,  we  have  undertaken  a  thorough  temperature  dependent 
study  of  solvation  and  rotational  dynamics  in  HEATN  using  C153  as  a  probe  and  compared  them  with 
results  previously  obtained  from  earlier  studies  on  imidazolium  and  phosphonium  ionic  liquids.  The 
solvation  time  in  HEATN  was  faster  than  the  isoviscous  imidazolium  and  phosphonium  ionic  liquids, 
which  is  probably  due  to  the  smaller  size  of  the  cationic  moiety.  The  solvation  times  in  HEATN 
showed  an  excellent  correlation  with  viscosity,  and  the  ionic  radii.  The  fluorescence  depolarization 
studies  showed  that  the  rotational  time  of  Cl 53  in  HEATN  at  different  temperatures  lies  within  the 
limits  of  hydrodynamic  predictions.  The  presence  of  the  amino  and  hydroxyethyl  groups  allow  HEATN 
to  participate  in  significant  hydrogen  bonding  interactions.  Surprisingly,  these  interactions  do  not  seem 
to  produce  a  substantial  difference  in  the  dynamics  of  the  HEATN  system  with  respect  to  imidazolium 
systems.  An  exception  is  given  in  Figure  6a,  where  HEATN  shows  a  linear  correlation  with  a 
significant  difference  in  slope  compared  to  the  imidazolium  and  phosphonium  systems.  Despite  the 
inherent  difference  of  having  the  additional  nitrogen  atom  in  the  ring,  no  significant  changes  in  the 
dynamics  of  this  system  were  observed.  Finally,  to  our  knowledge  this  the  first  study  of  solvation 
dynamics  in  a  triazolium  based  ionic  liquid. 
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Table  1.  Data  from  dimer  HEATN  structures  obtained  via  MP2. 


Dimer* 

Ereia 

(kcal/mol) 

OH-Ob 

o 

(A) 

NH-Ob 

o 

(A) 

Mulliken 
Charges 
(cation,  anion)e 

Geodesic 
Charges 
(cation,  anion)c 

& 

I 

0 

1.85 

1.90 

0.89,  -0.89 

0.82,  -0.82 

& 

0 

II 

0.6 

1.84 

1.96 

0.83, -0.83 

0.79,  -0.79 

£ 

0 

hi 

3.3 

1.87 

1.93 

0.80,  -0.80 

0.83,  -0.83 

IV 

3.6 

1.88 

1.91 

0.84,  -0.84 

0.80,  -0.80 

V 

4.0 

1.84 

1.93 

0.83, -0.83 

0.80,  -0.80 

w 

VI 

5.0 

1.86 

1.84 

0.88,  -0.88 

0.82,  -0.82 
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a  Erei  reports  the  energies  relative  to  dimer  “I”. 

b  OH-O  &  NH-0  refer  to  the  distance  between  the  OH/NH  group  on  the  cation  and  the  nearest  O  atom 
on  the  anion. 

L  Charges  on  each  cation  and  anion  are  reported  listing  the  cation  first  and  the  anion  second  for  both  the 
Mulliken  charges  and  the  geodesic  charges. 

*  MP2  optimized  dimer  structures  of  HEATN  using  the  6-31++G(d,p)  basis  arranged  in  order  of 
increasing  energy  with  “I”  being  the  lowest  energy  structure  found.  The  atoms  are  color  coded  such  that 
nitrogen  is  represented  by  blue,  oxygen  by  red,  carbon  by  black  and  hydrogen  by  white. 
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Table  2.  MP2  data  of  tetramer  HEATN  structures  (see  Figure  2) 


Tetramer 

Erda 

(kcal/mol) 

Ionb 

NH-Oc 

o 

(A) 

OH-Nc 

o 

(A) 

Mulliken 

Charge 

Geodesic 

Charge 

I 

0.0 

Cl 

N/A 

1.84 

0.88 

0.86 

A1 

-0.88 

-0.86 

C2 

N/A 

1.84 

0.88 

0.86 

A2 

-0.88 

-0.86 

II 

10.1 

Cl 

1.77 

1.70 

0.88 

0.86 

A1 

-0.93 

-0.74 

C2 

1.73 

N/A 

0.92 

0.78 

A2 

-0.88 

-0.91 

III 

18.8 

Cl 

1.93 

1.78 

0.95 

0.93 

A1 

-0.84 

-0.92 

C2 

1.88 

1.89 

0.73 

0.82 

A2 

-0.84 

-0.84 

IV 

27.8 

Cl 

1.89 

1.81 

0.76 

0.81 

A1 

-0.77 

-0.78 

C2 

1.86 

1.80 

0.70 

0.73 

A2 

-0.68 

-0.76 

V 

39.4 

Cl 

N/A 

1.73 

0.67 

0.95 

A1 

-0.79 

-0.90 

C2 

N/A 

1.75 

0.94 

0.83 

A2 

-0.80 

-0.88 

aErei  reports  the  energies  (kcal/mol)  relative  to  tetramer  “I” 

h  The  ions  are  named  as  follows:  the  Cl  ion  is  the  ion  that  is  perpendicular  to  the  plane  of  the  paper  & 
A1  is  the  anion  that  is  most  closely  associated  with  Cl  (see  Figure  2) 

c  OH-O  &  NH-0  refer  to  the  distance  between  the  OH/NH  group  on  the  cation  and  the  nearest  O  atom 
on  the  associated  anion. 
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Table  3.  Data  obtained  from  evaluating  the  MP2  tetramer  geometries  (see  Figure  2)  with  the  FM02- 
MP2  method  and  the  FM03-MP2  method. 


Tetramer 

Ere,a  (kcal/mol) 

Ionb 

FM02-MP2 

FM03-MP2 

MP2 

FM02- 

MP2C 

FM03- 

MP2d 

Mulliken 

Charge 

Geodesic 

Charge 

Mulliken 

Charge 

Geodesic 

Charge 

I 

0.00 

-0.30 

0.02 

Cl 

0.90 

0.93 

0.87 

0.83 

A1 

-0.90 

-0.93 

-0.87 

-0.83 

C2 

0.90 

0.94 

0.87 

0.83 

A2 

-0.90 

-0.93 

-0.87 

-0.83 

II 

10.13 

0.63 

0.02 

Cl 

0.83 

0.97 

0.88 

0.92 

A1 

-0.87 

-0.80 

-0.94 

-0.78 

C2 

0.90 

0.87 

0.93 

0.82 

A2 

-0.86 

-1.03 

-0.87 

-0.96 

III 

18.76 

0.55 

-0.09 

Cl 

0.91 

0.85 

0.93 

0.91 

A1 

-0.89 

-0.92 

-0.85 

-0.88 

C2 

0.91 

0.93 

0.81 

0.81 

A2 

-0.93 

-0.86 

-0.89 

-0.84 

IV 

27.75 

0.98 

-0.28 

Cl 

0.80 

0.80 

0.76 

0.75 

A1 

-0.89 

-0.78 

-0.79 

-0.73 

C2 

0.91 

0.74 

0.72 

0.68 

A2 

-0.83 

-0.77 

-0.69 

-0.71 

V 

39.43 

0.40 

0.13 

Cl 

0.81 

1.01 

0.67 

0.93 

A1 

-0.90 

-0.93 

-0.77 

-0.89 

C2 

0.95 

0.91 

0.91 

0.83 

A2 

-0.86 

-0.99 

-0.82 

-0.86 

aErei  reports  the  energies  (kcal/mol)  relative  to  tetramer  “I” 

b  The  ions  are  named  as  follows:  the  Cl  ion  is  the  ion  that  is  perpendicular  to  the  plane  of  the  paper  & 
A1  is  the  anion  that  is  most  closely  associated  with  Cl  (see  Figure  2) 

L  The  FM02-MP2  energy  is  relative  to  the  MP2  energy  for  the  same  tetramer 
d  The  FM03-MP2  energy  is  relative  to  the  MP2  energy  for  the  same  tetramer 
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Table  4.  Solvation  and  anisotropy  parameters  in  different  ionic  liquids. 


Ionic  Liquid 

T 

(K) 

(cP) 

R+b 

(A) 

R  1 

o 

(A) 

Anisotropyd 

Solvation 

ro 

<'Xrot> 

(ns) 

c  e 

'--rot 

/obs* 

<'Xsol\'> 

(ns) 

t+ 

(ns) 

t_ 

(ns) 

HEATN 

298 

427 

2.90 

2.09 

0.35  ±0.03 

21.1  ±0.62 

0.48 

0.36 

0.62 

4.7 

1.3 

308 

272 

0.36  ±0.10 

11.3  ±0.45 

0.42 

0.47 

0.22 

2.9 

0.70 

318 

133 

0.32  ±0.02 

5.7  ±0.38 

0.45 

0.45 

0.16 

1.4 

0.30 

328 

69 

0.38  ±0.02 

4.3  ±0.68 

0.67 

0.53 

0.10 

0.7 

0.18 

338 

40 

0.36  ±0.01 

2.6  ±0.08 

0.72 

0.60 

0.07 

0.4 

0.10 

BMIM  Cl1 

343 

334 

3.29 

2.09 

0.21 

0.59 

BIM  NTfV 

293 

90 

3.16 

3.39 

0.65 

0.20 

MIM  NTf2* 

293 

77 

2.69 

3.39 

0.64 

0.42 

CVIM  NTf2j 

342 

38 

4.05 

3.39 

0.62 

0.28 

P(C4)3C16  Brk 

339 

268 

4.96 

2.26 

0.29 

5.3 

P(C6)3Ci4  Clk 

331 

80 

5.18 

2.09 

0.33 

3.9 

a  Error  bars  for  viscosity  measurements  are  within  ±1-2%. 

van  der  Waals  volume  of  the  cations  (Vw)  are  calculated  from  the  Edward’s  atomic  increments  ,  and  the  radius  of  (R+)  are  obtained  using  the 
relation  R.  =  (3VW  /  4^) 1/3 

c  Radius  of  anions  (R_)  are  obtained  from  Bondi79  and  Jin  et  al.77 

d  The  error  bars  in  the  anisotropy  experiments  are  obtained  from  triplicate  measurements. 
e  Rotational  decoupling  constant  (Crot)  is  calculated  using  equation  4. 

Fractional  solvation,  /Gbs  are  calculated  at  50  ps  for  HEATN  samples  and  100  ps  for  other  ionic  liquids. 
g  Average  solvation  times,  <xsoiv>  are  obtained  by  the  fitting  C(t)  curves  with  two  exponentials. 

h  The  cation  (t+)  and  anion(t_)  diffusion  times  are  calculated  using  the  equation  t_  =  rjRl  /  f±kBT ,  where  the  parameters  are  defined  in  the  text, 
i  From  reference10 
j  From  reference68 
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k  From  reference12 
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Figure  captions 


Figure  1.  The  structures  of  (a)  solvatochromic  probe  Coumarin  153  (C153),  (b)  4-amino- 1,2,4-triazole, 
(c)  l-hydroxytheyl-4-amino-l,2,4-triazolium  nitrate  (HEATN),  (d)  1  -butyl-3 -methylimidazolium 
chloride  (BMIM+  Cl  ),  (e)  1-butylimidazolium  bis(trifluoromethylsulfonyl)imide  (BIM+  NTfi-),  (f)  1- 
methylimidazolium  bis(trifluoromethylsulfonyl)imide  (MIM+  NTf2~),  (g)  l-cetyl-3-vinylimidazolium 
bromide  (CVIM+  Br“),  (h)  tetradecyl  trihexyl  phosphonium  chloride,  and  (i)  hexadecyl  tributyl 
phosphonium  bromide  ( (C433C 1  nP+  Br~). 

Figure  2.  Tetramer  structures  of  HEATN  obtained  using  the  MP2  method  with  the  6-31++G(d,p)  basis. 
They  are  labeled  in  order  of  increasing  energy.  In  each  case  cation  1  is  the  cation  that  is  perpendicular  to 
the  plane  of  the  page  and  anion  1  is  the  closest  corresponding  anion.  Red  represents  oxygen,  blue 
represents  nitrogen,  black  represents  carbon  and  white  represents  hydrogen.  Additional  data  obtained 
from  these  calculations  can  be  seen  in  Table  3. 

Figure  3.  Representative  normalized  wavelength-resolved  fluorescence  decay  traces  of  Cl 53  at  three 
different  wavelengths  in  HEATN  at  55°C.  The  presence  of  the  decay  and  rise  components  as  a  function 
of  wavelength  is  a  key  signature  of  solvation  dynamics. 

Figure  4.  Time -resolved  fluorescence  emission  spectra  (TRES)  of  C153  in  HEATN  at  25°C  (top  panel) 
and  55°C  (bottom  panel)  at  50  ps  and  7.05  ns.  Corresponding  “steady- state”  and  “time-zero”  spectra 
are  also  included.  Note  that  the  7.05  ns  spectra  cross  the  “steady- state”  spectrum  at  25°C,  indicating 
that  these  spectra  do  not  adequately  model  the  completed  solvation  of  C153,  as  would  usually  be 
expected.  The  error  for  the  steady-state  and  time-zero  spectra  is  -  ±  100  cm  ',  while  for  the  other  timed 
points  it  was  ~  ±  200  cm4. 
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Figure  5.  Solvation  correlation  function  ( C(t ))  decay  curves  of  C153  in  HEATN  plotted  as  a  function  of 
temperature.  The  decays  were  fit  with  two  exponentials  and  the  average  solvation  time  as  well  as  the 
fractional  solvation  values  can  be  found  in  Table  4.  The  typical  error  for  each  data  point  is  <  0.05. 
Water  content  was  determined  to  be  <  0.3  wt.  %. 

Figure  6.  (a)  The  average  solvation  time  (<xs>)  of  Cl 53  is  plotted  as  a  function  of  viscosity  (q)  in 
different  ionic  liquids  based  namely  on  imidazolium,  phosphonium  and  triazolium  cations.  The 
solvation  times  of  C153  in  HEATN  show  a  reasonable  correlation  with  viscosity,  whereas  in 
imidazolium  and  phosphonium  ionic  lqiuids  studied  previously  poor  agreement  is  seen.  The  straight  line 
is  drawn  to  help  guide  the  eye.  (b)  Correlation  of  <xs>  in  the  ionic  liquids  with  q/T  and  cation  radius 
(R+).  The  proportionality  of  <xs>  varies  as  (i)/T)A  Rn  ,  where  R+  is  the  radius  of  the  cation,  and  A  and  B 

are  constants  determined  to  be  0.5  and  4.5  by  multiple  regression  analysis.  This  correlation  is  slightly 
improved  when  the  anion  size  (R .)  is  incorporated  (c).  These  plots  show  a  strong  dependence  of  the 
solvation  time  on  ion  size.  The  R"  and  standard  deviation  (SD)  are  obtained  from  the  fit  represented  by 
the  straight  line  in  ‘b’  and  ‘c.’ 

Figure  7.  Solvation  times  of  C153  in  HEATN  at  five  different  temperatures  versus  estimates  of  the 
cation  and  anion  diffusion  time  determined  using  the  equation  t±  =  T]Rl  /  f+kBT ,  where  the  parameters 
are  defined  in  the  text.  The  linear  fits  are  shown  as  solid  straight  lines,  whereas  the  dashed  line 
represents  the  1 : 1  correlation. 

Figure  8.  Correlation  of  the  average  rotaional  time  (<xrot>)  of  C153  in  HEATN  with  viscosity  (q/cP) 
and  temperature  (T/K).  The  solid  line  (red)  shows  the  best  fit  of  the  data  to  the  proportionality  <xrot>  = 

1.13xl0  9  (sKcP ’)  (q/T). 0  8  The  dashed  lines  are  the  average  rotational  times  predicted  from 
Distribution  A:  Approved  for  public  release;  distribution  unlimited 


OT 

hydrodynamic  models  assuming  an  ellipsoidal  shape  of  C153  using  “stick”  and  “slip”  boundary 
conditions.  The  experimentally  observed  rotational  times  fall  within  the  boundaries  of  hydrodynamic 
predictions  and  the  deviations  cab  be  estimated  by  the  rotational  coupling  factor  (Crot)  as  described  in 
the  text. 
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Figure  1. 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


Figure  2. 
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Figure  3. 
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Figure  4. 
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Figure  7. 


Figure  8. 
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